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Abstract 

Recent observations of the spectrum of J /if) produced in e + e~ collisions at the Y(4f>) resonance 
are in conflict with fixed-order calculations using the Non-Relativistic QCD (NRQCD) effective 
field theory. One problem is that leading order color-octet mechanisms predict an enhancement 
ON . of the cross section for J/ip with maximal energy that is not observed in the data. However, in 

this region of phase space large perturbative corrections (Sudakov logarithms) as well as enhanced 
nonperturbative effects are important. In this paper we use the newly developed Soft-Collinear 
Effective Theory (SCET) to systematically include these effects. We find that these corrections 
significantly broaden the color-octet contribution to the J ftp spectrum. Our calculation employs 
a one-stage renormalization group evolution rather than the two-stage evolution used in previous 
SCET calculations. We give a simple argument for why the two methods yield identical results to 
Q_i' lowest order in the SCET power counting. 
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I. INTRODUCTION 



Bound states of heavy quarks and antiquarks have been of great interest since the dis- 
covery of the J /if Q. In particular the production of quarkonium is an interesting probe of 
both perturbative and nonperturbative aspects of QCD dynamics. Production requires the 
creation of a heavy QQ pair with energy greater than 2mg, a scale at which the strong cou- 
pling constant is small enough that perturbation theory can be used. However, hadronization 
probes much smaller mass scales of order itiqv 2 , where v is the typical velocity of the quarks 
in the quarkonium. For J /if), rriQV 2 is numerically of order Aqcd so the production process 
is sensitive to nonperturbative physics as well. 

A systematic theoretical framework for handling the different scales characterizing both 
the decay andproduction of quarkonium is Non-Relativistic Quantum Chromodynamics 
(NRQCD) NRQCD solves important conceptual as well as phenomenological prob- 

lems in quarkonium theory. For instance, perturbative calculations of the inclusive decay 
rates for \c mesons in the color-singlet model suffer from nonfactorizable infrared diver- 
gences 0|. NRQCD provides a generalized factorization theorem that includes nonpertur- 
bative corrections to the color-singlet model, including color-octet decay mechanisms. All 
infrared divergences can be factored into nonperturbative matrix elements, so that infrared 
safe calculations of inclusive decay rates are possible [5j]. In addition, color-octet produc- 
tion mechanisms are critical for understanding the production of J /if) at large transverse 
momentum, p±, at the Fermilab Tevatron [6|. NRQCD has been applied to the production 
and decay of quarkonium in various experimental settings. However, there are still many 
challenging problems in quarkonium physics that remain to be solved Q. One important 
problem is the polarization of J /if) at the Tevatron. NRQCD predicts the J /if) should be- 
come transversely polarized as the p± of the J/ip becomes much larger than 2m c 0. The 
theoretical prediction is consistent with the experimental data at intermediate p±, but at 
the largest measured values of p± the J / if) is observed to be slightly longitudinally polarized. 
At these p±, discrepancies at the 3a level are seen in both prompt J /if} and if>' polarization 
measurements 0. 

New problems have arisen as a result of recent measurements of the spectra of J Jib 
produced at the T(AS) resonance in e + e~ collisions by the BaBar and Belle experiments [10|, 



111 ]. Leading order NRQCD calculations predict that for most of the range of allowed energies 
prompt J /if) production should be dominated by color-singlet production mechanisms, while 
color-octet contributions dominate when the J /if) energy is nearly maximal. Furthermore, 
as pointed out in Ref. [12], color-octet processes predict a dramatically different angular 
distribution for the J /if). Writing the differential cross section as 

S( PlP )(l + A(p^cos 2 9), (1) 



dpi/, d cos 9 

where p^ is the J /if) momentum and cos 9 is the angle of the J /if) with respect to the 
axis defined by the e + e~ beams, one finds the color-singlet mechanism gives A(p^) ~ 
except for large p^, where A{p^) becomes large and negative. On the other hand, color- 
octet production predicts A(p$) ~ 1. The significant enhancement of the cross section 
accompanied by the change in angular distribution were proposed as a distinctive signal of 
color-octet mechanisms in Ref. [12j . It was expected that these effects would be confined to 
J /if) whose momentum is within a few hundred MeV of the maximum allowed. 

Experimental results do not agree with these expectations. The cross section data as a 
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function of momentum is binned in intervals of 300 or 500 MeV, depending on the experi- 
ment, and the data does not exhibit any enhancement in the bins closest to the endpoint. 
However, the total cross section measured by the two experiments exceeds predictions based 
on the color-singlet model alone. The total prompt J/ip cross section, which includes feed- 
down from ip' and \c states but not from B decays, is measured to be a tot = 2.52±0.21±0.21 
pb by BaBar, while Belle measures = 1.47±0.10±0.13 pb. Estimates of the color-singlet 
contribution range from 0.4 — 0.9 pb [3, 3 IE 3|- Furthermore, A(j>$) is measured to be 



consistent with 1 (with large errors) for p^ > 2.6 GeV (Belle) and p^ > 3.5 GeV (BaBar). 
These aspects of the data suggest that there is a substantial color-octet contribution which is 
not confined to the very endpoint of the momentum spectrum but spread over a substantially 
broader range of momentum. 

In this paper, we will argue that higher order perturbative and nonperturbative correc- 
tions that are enhanced near the endpoint give rise to a broad color-octet contribution. 
The calculation depends on a nonperturbative function, called a shape function, which 
parametrizes the distribution of lightcone momentum of the J/ ip carried by the color-octet 
cc pair produced in the short- distance process. Since the shape function is nonperturbative, 
our calculation is not predictive. However, moments of the shape function are NRQCD 
operators whose size is constrained by the velocity scaling rules of NRQCD. Choosing a 
simple ansatz for the shape function whose moments are consistent with velocity scaling 
rules, we find that the combined perturbative and nonperturbative effects lead to substan- 
tial broadening of the color-octet spectrum in a manner that is consistent with data. Since 
the shape function that appears in this calculation also appears in other processes, it could 
be extracted from e + e~ data and used to make predictions for photoproduction and elec- 
troproduction once resummed calculations for these processes become available. Ref. [T3| 
presents a calculation of J/i/j photoproduction which includes a nonperturbative shape func- 
tion but not the resummation of large perturbative corrections. An important point of this 
paper is that nonperturbative and perturbative corrections combine to produce a spectrum 
that is broader than that obtained when only one of the corrections is included. Therefore 
it would be interesting to revisit the calculation of photoproduction using the methods of 
this paper. 

While the calculations of this paper show that the leading color-octet contribution is broad 
enough to be compatible with the observed distributions, other features of the e + e~ data 
remain puzzling. In particular, Belle reports a large cross section for J ftp produced along 
with open charm 



cx(e + e — > J/ipcc) 
a(e+e- -»• J/ipX) 



o-59±8:ii ±°- 12 



The predicted ratio from leading order color-singlet production mechanisms alone is about 
0-2 r! 3 an d a large color-octet contribution makes this ratio even smaller. In addition to 



the inclusive measurements, Belle reports a cross section for exclusive double charmonium 
production which exceeds previous theoretical estimates. Recent attempts to address the 
latter problem can be found in Ref. (l9| . 

The experiments also measure the polarization of the J/ip, which can provide impor- 
tant information about the production mechanism. Unfortunately the current experimental 
situation is unclear. Polarization is studied by measuring 

1 - 3?7l 

« = -j— , 

1 + Vl 
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FIG. 1: Leading order short distance amplitudes for e + e — » cc + g in QCD. 



where rjL is the fraction of J ftp which are longitudinally polarized. Both BaBar and Belle 
measure a ~ —0.5 for p^ < 3.5 GeV. However, for > 3.5 GeV BaBar measures a = 
—0.8 ± 0.09, corresponding to almost completely longitudinally polarized J/ip, while Belle 
measures a = —0.2 ± 0.2, which is consistent with no polarization. Neither measurement 
attempts to correct for feeddown effects on the polarization of the J/ip. Measurements 
of the polarization of directly produced J/tp can discriminate between various production 
mechanisms |l5j. The color-singlet process e + e~ — > J/ipcc produces J/ip with a rising from 
zero at p^ = to 1 at the kinematic endpoint for this process. The color-singlet process 
e + e~ — > J/ipgg prefers longitudinally polarized J/ip with a decreasing to almost —1 at 
the kinematic endpoint. Finally, the leading order color-octet diagrams produce J/ip with a 
between and -0.07, depending on the relative importance of S'-wave and P-wave production 
mechanisms. 

The NRQCD factorization formalism shows that the differential J/ip cross section can be 
written as 

da(e + e- -> J/V> + X) = da(e + e~ -> cc[n] + X) (O^) , (2) 

n 

where da is the inclusive cross section for producing a cc pair in a color and angular momen- 
tum state labeled by [n] = 2S+l L <y j\ In this notation, the spectroscopic notation for angular 
momentum quantum numbers is standard and i — 1(8) for color-singlet (-octet) production 
matrix elements. The short-distance coefficients are calculable in a perturbation series in 
a s . The long-distance matrix elements (On ) are vacuum matrix elements of four-fermion 
operators in NRQCD 0. These matrix elements scale as some power of the relative velocity 
v <C 1 of the c and c quarks as given by the NRQCD power-counting rules. 

At lowest order in v the only term in Eq. (J2J) is the color-singlet contribution, [n] = 3 S^\ 
which scales as v 3 . The coefficient for this contribution starts at 0(a 2 s ) Away from the 
kinematic endpoint E max — (s + M|)/(2\/s), where s is the center-of-mass energy squared, 
color-octet contributions also start at 0(a 2 s ). Since the color-octet contributions are sup- 
pressed by f 4 0.1 relative to the leading color-singlet contribution, they are negligible 
throughout most of the allowed phase-space at leading order in perturbation theory. How- 
ever, as pointed out in Ref. 12||, there is an 0(a s ) contribution to color-octet production 
from the diagrams shown in Fig. ^ In this process the cc are in either a 1 S$ S) or 3 Pj 8) 
configuration. The resulting cross sections are proportional to a 5— function which forces 
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the energy of the cc pair to be maximal l3| : 

az a cos a s z m c 

da[ 3 p( 8) ] _ _ 8vr 2 a 2 a s e 2 (3 + 6r + 7r 2 ) + (3 - lOr + 7r 2 ) cos 2 
dzd cos 9 s 2 m? c 2(1 — r) 

where r = 4m 2 /s, and z = E c5 /E^ ax with -E^ aa: = a/s(1 + r)/2. Here we have used 
(Og( 3 Pj)) = (2 J + l)(Cg ( 3 P )) and combined the short- distance cross sections for the 
three P-wave contributions. For comparison, in this region of phase space the color-singlet 
contribution approaches a constant 

.. da[ 3 sP] 647ra 2 a 2 e 2 i .fl + r 2 \ 

hm - , ' = -— r- — (1 + r) cos 2 6 . (4) 

z^idzdcos6 27s 2 m c v ; ll-r J w 

We see that the v 4 suppression of the color-octet contribution is compensated by one power 
of a s . In addition, the numerical coefficient of the color-octet contributions is bigger by a 
factor of ~ 47r which is due to the difference between two- and three-body phase space. Since 
the color-octet contribution is proportional to a singular distribution, sensible comparison 
with data requires that we integrate the cross section over a region in z of finite size. It is 
important to know what range of z must be integrated over in order for the leading order 



calculation to be reliable. Ref. 21] has identified a class of nonperturbative corrections which 
are suppressed by powers of v 2 but are enhanced by factors of 1/(1 — z), so we must integrate 
from a lower limit no larger than z m i n ~ 1 — v 2 . For z min ~ 1 — v 2 , the integrated color- 
singlet cross section gets an additional v 2 suppression relative to the color-octet contribution 
so the ratio of the integrated color-octet to integrated color-singlet is ~ 4iTV 2 /a s which is 
roughly a factor of 10. 

Note that because v 2 ~ 0.3 for J/ip, the cross section must be integrated over a substantial 
region. For production at the T(4S), E max = 5.74 GeV, corresponding to p m ax = 4.84 
GeV, so we must integrate above E w 4.0 GeV, which corresponds to ~ 2.5 GeV. 



These numbers show that the spectrum reported in Refs. jlOl . Illj . which plot the data as 
a function of p^, cannot be directly compared to the leading order NRQCD calculation. 
Nonperturbative effects smear the color-octet contribution over nearly half the range of 
allowed momenta, so the color-octet contribution is not confined to the last momentum 
bin. Higher order perturbative corrections are also important. The next-to- leading order 
radiative correction to color-octet production has a contribution proportional to a s log(l — 
z)/(l — z), so perturbative corrections require a resummation when 1 — z ~ a s . Since 
a s ~ v 2 , similar conclusions regarding the reliability of the leading order calculation arise 
from purely perturbative considerations. 

Large perturbative and nonperturbative corrections signal that the NRQCD factorization 
theorem of Eq. (J2J) is no longer valid. The problem is that near the endpoint the J/ip is 
recoiling against a gluon jet with energy of order M T but mass of order M T a/ Aqcd / . 
These degrees of freedom are integrated out of NRQCD, but when we probe the endpoint 
region they must be kept as explicit degrees of freedom. The effective theory which correctly 
describes this kinematic regime is a combination of NR QCD for the heavy degrees of freedom, 
and the Soft-Collinear Effective Theory (SCET) (22I I23LI24I lij^ for the light energetic degrees 
of freedom. The factorization theorem derived using NRQCD and SCET will include a 
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nonperturbative distribution that incorporates the v 2n /(l — z) n nonperturbative corrections 
to all orders. In addition, renormalization group equations of SCET can be used to resum 
large perturbative corrections. A similar treatment of nonperturbative and perturbative 
endpoint corrections to the color-octet contributions in the inclusive decay T — > X + 7 can 
be found in Ref. 

II. FACTORIZATION 

In this section, we present the factorization theorem for J/ip production near the kine- 
matic endpoint. We begin by studying the kinematics of the process to determine when 
NRQCD breaks down and SCET must be used instead. In the e + e~ center-of-mass (COM) 
frame, the cc pair has momentum p^ E = Mv M + i^, where M = 2m c , ^ is the residual 
momentum of the cc pair and the four- velocity of the J/ip is 



Here is the J/ip mass and x = (E^ + p^)/\fs. The residual momentum arises because 
the cc pair is produced along with ultrasoft (usoft) quanta which carry O(Aqcd) momentum 
in the rest frame of the J/ip. The lightlike vectors are n M = (1, 0, 0, 1) and n M = (1, 0, 0, — 1), 
where the J /if) is moving in the z direction in the COM frame. It is sometimes helpful 
to write I 11 in terms of the boost that takes one from the J/ip rest frame to the COM 
frame: £ M = A M ^. The components of are O(Aqcd) while the components of £ M scale 
as: n ■ £ ~ M^Aqcd/(^V^)> n ' & ~ ^a/sAqcd/-^ and £± ~ Aq CD . The momentum of the 
virtual photon is = ■\/s/2(n fl + n M ) and the gluon jet has momentum 



where x = xM/M^. Away from the endpoint region (x 1) the recoiling gluon jet momen- 
tum and invariant mass are both of order s/s. Therefore, the jet can be integrated out of 
the theory and one obtains the NRQCD factorization formula in Eq. (J2J). In this region £ M is 
negligible compared to the large components of Px and p% and can be set to zero. Therefore 
the cross section is not sensitive to motion of the cc within the J/if>. This is evident from 
the NRQCD factorization formula which depends on a single nonperturbative parameter, 
(O^Sf } )) or (C( 3 P (8) )). When 1 - x ~ A QCD /M ~ v 2 , y/s(l - x) is the same size as 
n ■ £, so the residual momentum cannot be neglected. In the endpoint region, the NRQCD 
factorization formula breaks down. A new factorization theorem is needed which includes 
a distribution function that parametrizes the nonperturbative motion of the cc pair within 
the jet containing the J/if>. For 1 — x = 1/3, p^p = 2.8 GeV, so the new factorization theo- 
rem is relevant for a significant part of the measured p^ spectrum. The failure of NRQCD 
factorization can also be understood by considering the gluon jet which was integrated out. 
When 1 — x ~ Aqcd/^ the jet is no longer highly virtual: 



Since m x /E x ~ Aqcb/M <C 1, the gluon jet is composed of energetic particles with small 
invariant mass that must be included explicitly in the effective theory. 




(5) 





m 2 x ~ Aqcds/M . 



(7) 
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SCET has collinear degrees of freedom whose momentum scales as n -p ~ Q, n -p ~ A 2 Q, 
and p L ~ XQ. For the process under consideration, Q is of order y^, while A ~ \Jl — x ~ 
(Aqcd/M) 1 / 2 . SCET also has soft de grees of freedom whose momentum scales as A and usoft 
degrees of freedom whose momentum scales as A 2 . Heavy quark fields in SCET are the same 
as in NRQCD when considering quarkonium. Since a field redefinition removes the large 
0(m c ,m c v) part of their momentum, derivatives on these fields scale as Aqcd and therefore 
they are usoft fields in SCET. Thus the endpoint region of J/ip production is mediated by 
SCET operators involving usoft (heavy) quark fields and collinear gluon fields. Soft fields 
do not enter to the order we are working, so are neglected. 

To match onto SCET, matrix elements in QCD are evaluated at the scale Q and expanded 
in powers of A. Each order in the A expansion is reproduced in the effective theory by a 
product of Wilson coefficients (which depend only on the large scale Q) and SCET operators. 
At each order, one must include all SCET operators which can contribute to the process 
under consideration, subject to the restriction that the operators respect the symmetries 
of the effective theory. In SCET, gauge transformations can be classified by their scaling 
with A just like fields j25[. For e + e~ — > J/ip + X, the relevant operators must be invariant 
under collinear and usoft gauge transformations j25[. SCET is not Lorentz invariant, but the 
Lorentz invariance of QCD is realized in the form of constraints on the possible operators 
which appear in the Lagrangian and in matching calculations. These constraints are called 
reparametrization invariance (RPI) Usoft and collinear gauge invariance and RPI 

uniquely fix the form of the lowest order operators contributing to e + e _ — > J/ip + X. 

In the collinear sector of SCET there is a collinear fermion field £ njP , a collinear gluon 
field A% „, and a collinear Wilson line 



W n (x)= expi-g s ^n- A n:q (x)j . (8) 

L perms ^ ' - 

The subscripts on the collinear fields are the lightcone direction n^, and the large compo- 
nents of the lightcone momentum (n ■ q,q±). The operator V 1 projects out the momentum 
label . For example n ■ V^ n ^ = V£, n ,p — n • p £ UiP . Likewise in the usoft sector there is a 



usoft fermion field q us , a usoft gluon field A^ s1 and a usoft Wilson line Y . Using the trans- 
formations for each of these fields under collinear and usoft gauge transformations given in 
Ref. [25], we can build invariant operators. The collinear-gauge invariant field strength is 

GT = --Wi[iV% + g s A»iVZ + g^ n jW, (9) 
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s 



where 

WZ = —P + Vt + —in-D, (10) 

and iD^ = id^ 1 + g s A^ s is the usoft covariant derivative. RPI invariance requires the label 
operators and the usoft covariant derivatives, which scale differently with A, to appear in 
the linear combination appearing in iT>£. The leading piece of is order A and can be 
written as n u G v ^ = i[P, Bf], where 

B» ± = -W^V»+g s (A», q )±W. (11) 
9s 

The subscript _L on B^ indicates that /i must be a perpendicular direction. The leading 
operator for e + e~ — ► ccg must have a (A% )j_ field operator to create the collinear gluon in 
the final state. B*t is the collinear gauge invariant generalization of this field. 
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+ crossed diagram 



FIG. 2: Matching the production amplitude for e + e 
are represented by a spring with a line through it. 



cc + g in QCD and SCET. Collinear gluons 



The 0(A) operator which creates a heavy quark and antiquark in a 1 S state with a 
collinear gluon is 

^ p T^ s °\-n-vV^)BZx-p. (12) 

This is collinear gauge invariant because B± is invariant and the heavy quark fields do not 
transform under collinear gauge transformation. It is also usoft gauge invariant since 



X- 



V w x-p, ^ P -> V w ipp, B\ 



VusBlVl, 



(13) 



under usoft gauge transformations. Type-Ill reparametrization invariance |27[ requires P 
to appear multiplied by n ■ v j28|]. The operator is not invariant under type-I and -II 
transformations but can be made invariant by adding A suppressed operators which are not 
needed at the order we are working. The factor of — P in Eq. ()12|) gives the large component 
of the momentum of the jet at the endpoint. This is the part of n ■ px that survives as 
A-qcd/M and 1 — x — > 0, so that — n • vP = s(l — r)/M. At leading order, the function 

ri 8 ^ is determined by requiring the SCET matrix element of Eq. (|12|) to agree with the 
lowest order QCD diagrams for e + e~ — *> ccg, shown in Fig. 



M 



-a/i ' 



where e 



a/i 



rv 



p n f3 /2. The leading SCET color-octet 3 Pj operator is 



^A 3 f/\-n ■ vP; rfBZA ■ p CT A • <r 5 X - P , 
where p = p/M and the leading order matching coefficient is 



(14) 



(15) 



<8, 3 Pj) 



M 



9o 



1 + r 



where g 



~9aS 



1 - 



9^5 



-(1 - r)n s n, 
r 



-(1 + r)n^n c 



9 a ^ n s n ^ + n & n <T, 



(16) 
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It is possible to show that the rate e + e — > J /if) + X in the endpoint region can be 
factored into a hard coefficient o"o, a collinear jet function J, and a usoft function S, as 



where z = E^,/E^ ax , and 



P[r, z] 



a/(1 + r) 2 z 2 - 4r 



(17) 



(18) 



is a phase space factor. Note that P[r, 1] = 1. The derivation is given in Appendix El For 
convenience we have integrated over 9. The 9 dependence of the resummed cross sections 
are the same as in Eq. (jSJ). 

The jet function in Eq. (JTTj) is independent of the state of the cc pair in the J /if), and is 
defined as 



J(n ■ pn ■ k + p 2 ± ) 



(19) 



s(l + r) 
Air 



Im 



d A y e ik ' y (0|T (Tr 



(y) 



Tr 



'rMi(o) ) |o> 



where the factor s(l + r) is chosen to provide a convenient normalization for the process we 
are interested in. J(n • pn ■ k + p^_) is a function of the label momenta n ■ p and p 2 ± , and the 
usoft momentum n ■ k, which follows from the collinear Lagrangian containing only the n ■ d 
derivative j2^|. Furthermore RPI requires that the jet function depends on the combination 
n ■ pn ■ k + p\. Here we parameterize the jet momentum so that n ■ p = —n ■ px and p± = 0. 
The jet function is perturbatively calculable since it is determined by physics at the scale 
a/ sAqcd/^c ^> Aqcd- The 0(a s ) calculation of the jet function, as well as the derivation 
of the leading order renormalization group equations, is given in Appendix IBl 

The shape functions, defined in terms of usoft fields, depend on the state of the cc pair, 
so there are different shape functions for the and 3 -Pj 8 ' ) contributions. They are 



5(8.^0)^ 



IxlpT^p a\a^ 5(i + - in ■ D)^ p T B X - P \0) 
4m c (Ot(iS )) 

\ X l p Hg ■ oT B ) iP P ala^5(£ + - in ■ D)^ p (-jD ■ aT B ) x _ p |0) 

4m c <Oj(3p )> 



(20) 



(21) 



where £ + =n-£. Note that the combination d£ + S^(£ + ) is boost invariant, so d£ + S^(£ + ) = 
d£ + S^(£ + ), where £ + is the residual momentum in the J /if) rest frame. It is useful to go 
to the J /if) rest frame when evaluating the shape functions, since in that frame derivatives 
acting on the heavy quark fields scale as m c v 2 . The normalizations of these functions are 
defined so that J d£ + S^ 8, +lLj \£ + ) = 1. The jet function normalization is chosen so that 
the hard coefficients correspond to the lowest order total cross section for the color-octet 
contributions: 



(8,^0) 



32tt 



3s 2 - y 



— r 



(8, 3 Po) 




^(3 



32ir 2 a s a 2 e 2 c (Of ( 3 P )) 3 + 2r + 7r 2 



3s 2 



(22) 
(23) 
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To see how the factorization theorem in Eq. (|T7|) reduces to the leading order cross section 
in the appropriate limit, we use the result in Appendix |B] for the lowest order jet function, 

(24) 



where £ is defined by 



s(l + r) 



z + 









1 — r 


{-- 




1 + r 




M J 



= s(l + r)(£-z), 

and A = - M. Then we combine Eqs. (|I7jl . (1271 and 

is defined as 

S [n] (0 = M^-S [n] (£ + ) 
so that / rf£5 M (0 = 1, to obtain 



4m e (Of( 1 5 )> 



z — 



(25) 

, where as a function of £ 

(26) 



1 — r \ in ■ D — A 



1 + r 



M 



^ P T B X-p|0> 
(2?) 

for the 1 £q contribution, for example. Here, n ■ D = M^/(x-s/s) n ■ D and scales like m c v 2 
or Aq CD . When 1 — z ^> (m • D — A)/M, we can drop the derivatives inside the 5-function 
and pull the 5-function outside the matrix element. The final result is da( 8 ' 2S+lLj * > /dz = 



a, 



(8,^+i Lj ) 



6(1 — z), which agrees with Eq. (JHJ). 



III. RESUMMING SUDAKOV LOGARITHMS 



The leading order color-octet contribution is proportional to 6(1 — z). The next-to- leading 
order radiative corrections have contributions of the form a s ln(l — z)/(l — z). Clearly when 
z ~ 1 — a s these corrections are large and must be resummed. This can be accomplished 
in a straightforward manner by using the renormalization group equations of SCET. The 
resummation is most easily carried out by taking moments with respect to z, then the large 
corrections as z — ► 1 become large logs of N in the expression for the Nth. moment. In this 
section, the resummation for the color-octet 1 Sq production mechanism will be described 
explicitly. The moments of the cross section in Eq. (JT7jl are 

J\ zz n^_!1 = af Sa \fi)j\zz N j\tS^ 1So \Z,fi)J(s(l + r)(Z-z),fi) (28) 

= 4 8,15 °V) f\ f dzz N S^ ls °\{,»)J(s(l + r)(t-z),ri 
Jo Jo 

= <4 8,15 °V) / d££ N+1 [ duu N S( 8 > 1So \t,riJ(s(l + r)Z(l-u),n). 
Jo Jo 

In the last line of Eq. (J28)) . we have made the substitution z — «£. Since the large logs come 
from the region £, z ~ 1, P[r, z] can be replaced with P[r, 1] = 1 and the factor of £ in the 
argument of the jet function can be set equal to 1. Then the moments factorize: 

ap 1S ^ = *f S °\riSp 1S °\riJM, (29) 
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where 



a 



N 



s 



(8/50) 



N 



j n (ij) 



1 r/rT^ 8 ' 150 ) 
1 



dwu^ J(s(l+r)(l -«),//). 



(30) 



Note that we are only interested in the large N moments, so have used SV+i = Sn+0(1/N). 

To resum logarithms we must find the renormalization group equations for the three 
terms on the right hand side of Eq. (ffijj) . o"q 8 ' S °\l^) is proportional to the square of the 
matching coefficient for the lowest order 1 S'q 8 ^ operator Eq. (fT2"j) . so its anomalous dimension 
is simply twice the anomalous dimension of that operator, which is calculated in Ref. (2^ 



d 

d\x 



(8^0)/ X 



2C A a s 



71 



log 



MfL 
s(l — 7 



2a, 



77 



17 r n f 



(8^0) 



(/i) , (31) 



where fj 2 = 47re~ 7 /i 2 . erg 8 ' P °\^) has the same anomalous dimension. The scale appearing 
in the logarithm is — n ■ vV = s(l — r)/M. 

The renormalization group equations for the moments of the jet function are calculated 
in Appendix El The result is 



djj, 



2C A a s 



7T 



log 



fi 2 



N 



s(l +r)N 



+ 



2a, /ll 



7T 
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6 



J N (/i) 



(32) 



where N = e 
in the definition of Jn is s(l + r) 
anomalous dimension in Eq. (1321) . 



From Eq. (jSUj) it is clear that the only dimensional quantity appearing 



Since the moments a 



(S/So) 



so this is the scale that appears in the logarithm of the 



are physical and therefore ix independent, the renormalization 



group equation for Sn immediately follows: 



d\i 



2C A a s 



71 



log 



p, l-r N 
Ml + riVo 



+ 



a,C/ 



77 



S N (/i) 



(33) 



The anomalous dimension of Sn does not depend on the angular momentum quantum 
number of the cc state, so the resummation is identical for both 1 5'o and 3 Pq production 
mechanisms. Note that the scale M(l + r)/(l — r) appearing in Eq. (}3*3*|) also appears 
dividing in ■ D in Eq. (}2Tj) . 

We are able to obtain the anomalous dimension for Sn because we can calculate the 
renormalization group equations for the hard coefficient and jet function, and because the 
moments of the cross section are n independent. Direct evaluation of the anomalous dimen- 
sion for the shape function from the definition in Eq. (J2(Jj) is complicated by the projection 
operator a^a^. In quarkonium decays there is no such projection operator, and the deriva- 
tion of the evolution equation for the decay shape function is straightforward \2oL |26( . The 
result is remarkably similar to Eq. (J3*3*j) . The only difference is that the scale appearing in 
the logarithm is 2rriQ in decay as opposed M(l + r)/(l — r) above. Thus the coefficient of the 
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logarithm, often referred to as the cusp anomalous dimension, and the second term in the 
square brackets of Eq. (}3"3"j) is the same in production and decay. This holds to all orders a s 
because both processes obey a similar factorization theorem. Specifically the jet functions 
are identical in the two processes, and the hard functions obey the same renormalization 
group equation since the same SCET operator mediates production and decay. The only 
differences between the renormalization group equations for the two processes are the scales 
appearing in logarithms. 

Defining \xn = (s/M)(l — r) and yo = r(l+r) /(l — r) 2 (No/N), we see that the logarithms 
in the hard, jet and usoft functions are minimized at the scales fiHy/Vo and HhUo respec- 
tively. Large logarithms of N are resummed by evolving the jet and usoft functions to their 
respective scales. The evolution can also be done in one step by defining separate renor- 
malization scales for collinear and usoft loops. Loops whose momenta scale like (1, A 2 , A) 
come with a factor of \x^ D and loops whose momenta scales like (A 2 , A 2 , A 2 ) come with a 
factor Hu~ D - This idea is similar to the velocity renormalization group in NRQCD 0. The 
renormalization group equations for Jn and Sn take the form 



Hc-i — J n = 1j(^c)Jn, 

Hu-. — S N = js(n u )S N . 
dfi u 



(34) 



Factorization of usoft and collinear degrees of freedom guarantees that 7j is a function of 
fi c only and that 75 is a function of fi u only. The scales are however correlated, so that 
A*c = fJ'Hy/y and \i u = \ihV- Evolving the variable y from 1 to y simultaneously resums 
large logs in both J N and Sjy. 

Defining 1^ = JnSn, the evolution equation for T N as function of y is 



dy 



1 



7j(/W^) + 75 V (/^))r 



N 1 



N 



(35) 



This equation is easily integrated to obtain the following expression for the resummed mo- 
ments: 

(36) 



ff W8o) = ff f. 1 *) jS g. 1 *)( MflWj ) e «o B (N)p 1 (x)+»( X ) j 



where 

9i{x) 



2Ca 



[(1 - 2 X ) log(l - 2 X ) - 2(1 - x) log(l - x)\ 



(37) 



sr 



adj 



$1 



[- log(l - 2 X ) + 2 log(l - x)\ ~ log(l - X) 

1 



PI 
2C A 



log(l - 2 X ) - 21og(l - x) + ^log 2 (l - 2 X ) - log 2 (l - x) 



AC a ( 1 + r 
log(l -2 X )--^ log r- -N ) [log(l - 2 X ) - log(l - X )\ , 



Po " /W Po ~° V (1 - r) 2 

X = log(JV)a s (// H )A)/47r, Tf = C A [C A (Q7/3Q-n 2 /12) -5^/18], (3 = (UC A -2n f )/3, and 
A = {34C 2 A -lOC A n f -6C F nf)/3. Tf istheO(a 2 ) piece of the cusp anomalous dimension, 
which was taken from Ref. j29(. 
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The expression in Eq. ()36|) gives the resummed expression for the moments of the dif- 
ferential cross section to next-to-leading logarithmic order. To obtain the differential cross 
section, the inverse-Mellin transform of Eq. ()36|) must be taken. Using the results of Ref . [30] , 
we find: 

= -£fp[r,z]*W*) S V8o ){(:) (38) 
z—\e(i ~ N exp[lg 1 [a s (3 l/(47r)]+ g 2 {a s f3 l/(4n)]] 



dz{^ ' T[l - gi [a s [3 l/(4ir)} - a s (3ol/(^)g[[a s /3ol/^7r)]] 

where I ~ — log(£ — z), a s = a s (fin), and the shape function contains no large logarithms. 1 
To obtain the octet 3 P contribution let a^' 1 ^ -> a( 8 - 3p °) and S^ 1 ^ S^ 8 ' 3 ^. 

In this paper, the jet function and soft functions are first factorized and then run down 
to their appropriate scales. By linking the scales through the introduction of the variable 
y, large logarithms have been resummed in a single step using Eq. (J3*5j) . This approach to 
renormalization group evolution, which we will refer to as one-stage running, is identical to 
the method employed in Ref. 32] for the resummation of b — > sj. Early applications of 
SCET 113] adopt a two-stage running approach. First, operators in SCET are evolved from 
the hard scale, ^h, to the intermediate scale, which for the process under consideration is 
HHy/Vo- Then collinear degrees of freedom are integrated out, leaving only usoft degrees of 
freedom. Nonlocal operators in the usoft effective theory are then further evolved to the scale 
/ihVq- This approach to resummation in the process b — > sj has been shown to yield 
identical results at the next-to- leading log accu racy as one-stage running j32|. Applications 
of SCET to quarkonium decay in Refs. jiy, 3^, 34( also use a two-stage approach. 



To see why the two methods yield equivalent results (to leading order in A) we first note 
that_ the two-stage running approach of Ref. [I^f leads to the following evolution equations 
for r^v = JnSn'. 

V-t-Tn = (7j +7g)fjv VHy/Vo < A* < I*h (39) 



djj, 



-2 lH T 



N 



d vi 



ti dfjF N = 1sTn ^ Hy ° - ^ - VHy/yo- 

The first stage corresponds to running the coefficient of the SCET operator from the scale 
jin down to the scale (AHy/Vo- Note the evolution in this stage does not depend on the 
moment N. The second corresponds to the evolution of the usoft function in the purely 
usoft theory down to the scale HhDo- The difference in the two-stage running in Eq. (J39"J) 
and the one stage running in Eq. (j35j) can be visualized geometrically as integrating Eq. (|34jl 
along two different paths in the (/i u ,/i c ) plane j35[. The evolution in Eq. (|3*9*|) corresponds 
to the path 

(/i«,/i c ) = fm(y,y) Vm<y<^ (40) 

= m(y, Vm) y <y<y/]fy, 



1 The Landau pole in Eq. should be dealt with in the same fashion as in B decays |31|. 
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while integrating Eq. corresponds to the path 



(jmhUc) = f*H{y,Vy) 2/0 < 2/ < i • (41) 

Since the paths begin and end on the same point the difference can be expressed as an 
integral over a closed loop in the (/i u , /i c ) plane, and will be zero if the anomalous dimension 
vector, 7 = (jgjjj) has a vanishing curl: V x 7 = where V = (fi u d/dfi u , fi c d/dfj, c ), or 
equivalently, 

^ = ^- < 42 > 

At leading order in A, this equation is trivially satisfied due to factorization, so the two- 
stage and one-stage evolution will give identical results. This may not be true for operators 
appearing at higher orders in A since the factorization only holds to lowest order in A. 



IV. PHENOMENOLOGY 

Before we can investigate the phenomenological consequences of our analysis we must 
determine the shape function. Unfortunately not much is known about the octet shape 
functions in J/ip production. They also arise in photo- and electroproduction j2l|, so once 
the resummed calculations of these processes are available the universality of the shape 
functions could be tested. In this paper we use a model of the shape function to fit the 
available data so our calculation is not predictive. However, the shape function is not 
completely arbitrary because the moments of the shape function are NRQCD operators 
whose sizes are constrained by NRQCD power counting rules. For example, taking moments 
of Eq. ((201) gives: 

sf So) = j di + (i + ) N s^ 1So \i + ) 

= (OlxtpT^ajfrMm ■ D) N ^ p T A X- P \0) 
4m c (0j(%)> 

Each derivative is of order m c v 2 ~ Aqcd ; so the N th moment is 0(Aq CD ). 

For our model of the shape function we adopt a modified version of a model used in the 
decay of B mesons |36j], 



1 n ab P + 

f(£+) = IT py(* - lr^e-^" 1 ) , x = "-, (44) 

where a and b are adjustable parameters and A = — M . Note i + is the residual 
momentum of the cc pair in the rest frame of the J/ip. The first three moments of Eq. (J44|) 
are 

m = dt f(i + ) = 1, m x = / dt £+ f(£ + ) = A(b + 1), 

J A J A 

m 2 = dt (£ + ) 2 f(£ + ) = A 2 Q + (b + lA . (45) 



A 
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Since A ~ C(Aqcd), any choice with a ~ b ~ Oil) gives the desired scaling for the 
moments. For the sake of simplicity, we assume that the parameters in the shape function 
model are the same for both of the color-octet contributions. The data from both BaBar 
and Belle include feeddown from ip' and Xcj so these contributions must be taken into 
account. Because the Xc are P-wave states, the cross section for their production is v 2 
suppressed relative to J/tp and the feeddown to J /if) is further suppressed by relatively small 
branching fractions. For this reason, we neglect feeddown from \c m our analysis. However, 
feeddown from ib' states is not suppressed and must be included. The shape function for 
if}' can be different from that of J/ib but for simplicity we will assume the same form for 
both shape functions. Feeddown will affect the overall normalization of the color-octet 
cross section which is given by the linear combination (O^ ( 1 5'o)) + 3.8(Of ( 3 Po))/m 2 , where 
(Ot{ 2S+1 Lj)) = (C 8 J/ ^( 25+1 Lj)) + Br(^ -> J/if)){Of { 2S+1 Lj)). Likewise, the normalization 
of the color-singlet contribution is given by 

(Of{ 3 S 1 )) = (O^CSJ) +Br(V' -> J/if)){Of\ z S 1 )) = 1.45 GeV 3 . (46) 

In Fig. El we show as the solid line the shape function in Eq. (|44|) convoluted with the 
perturbative resummed result, normalized to a . The dashed line is a plot of the shape 
function alone, and the dotted line includes only the perturbative resummation. Here we 
use m c = 1.4 GeV, the one loop coupling constant is evaluated at the hard scale \in = 
(s/M)(l — r), with Aqcd = 190 MeV, and a = 1 and 6 = 2. The value of the first and second 
moments of the shape function for this choice of parameters are 890 MeV and (985 MeV) 2 
respectively. We have chosen the shape function parameters to give a reasonable description 
of the p^ distributions measured by BaBar and Belle. Since m c v 2 ~ 500 MeV the moments 
are consistent with the velocity scaling rules. Fig. El gives a picture of the effects of both the 
perturbative resummation and the shape function. The perturbative resummation gives a 
result that is highly peaked in the endpoint region. The shape function is also peaked close 
to the endpoint, though at a lower value, and is broader. It is interesting to note that the 
convoluted result is broader yet, and its peak is shifted to the left of both the perturbative 
resummed peak and the shape function peak. 

In order to make a consistent comparison of theory to data one needs to treat the endpoint 
of the color-singlet contribution in SCET and NRQCD. We leave this for future work, and 
will use the leading order NRQCD calculation of the color-singlet contribution over the 
full range of momenta. In Fig. 0] we show the sum of the color-octet and color-singlet 
contributions as the upper line, and the color-singlet contribution only as the lower line. 
The color-octet matrix elements which set the normalization are chosen to be <C^(%)) = 



(^( 3 P ))/m c 2 = 1.3 x lO" 1 GeV 3 . This is plotted against the BaBar data [ll(. In Fig. 
we show the same comparison to the Belle data [10]. The shape function parameters are 
the same, only the color-octet matrix elements are smaller: (Of ( 1 5'o)) = (Of ( 3 Po))/m 2 = 
6.6 x 10~ 2 GeV 3 . These values of the color-octet matrix elements are consistent with data 
from photo- and hadroproduction j3?J . Figs. 0] and El concisely present the central point 
of this work: when perturbative logarithms are resummed and the non-perturbative shape 
function is included the color-octet contribution becomes a broad distribution as a function 
of the J /if) momentum. Furthermore, for a reasonable choice of parameters including the 
color-octet contribution gives a fairly good description of the p^ distribution. The color- 
singlet contribution alone cannot describe the data on J/ if) production at BaBar and Belle. 

The effect of the resummation on the angular distibution of the J Jib is shown in Fig. |U] 
where we plot A(p^). The circles with error bars are data from Belle [lOj while the triangles 
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FIG. 3: The color-octet contribution at the endpoint normalized to ctq. The solid line is the pertur- 
bative resummation convoluted with the shape function. The dotted line is perturbative resummation 
only, and the dashed line is no resummation, but including a shape function. 

with error bars are data from BaBar . Note that the central value of one data point lies 
outside the physical range, —1 < A(p^) < 1. The solid line is the color-singlet prediction 
and the shaded band includes the resummed color-octet contribution. The color-octet x Sq 
production mechanism gives A = 1 while the 3 Pj mechanisms gives A « 0.7. The band 
shown in Fig. |H1 is obtained by varying the two matrix elements (Of^So)) and (C^( 3 P )) 
subject to the constraint that both matrix elements are positive and that 

(O*( l S )) + 3.8(Ot( 3 P ))/m 2 c = 3.2 x 10" 2 GeV 3 , (47) 

which is the normalization of the curve shown in Fig. |S| If the normalization in Fig. 0] 
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FIG. 4: The sum of the color-octet and color-singlet contributions are plotted as the upper line. The 
lower line is the color-singlet contribution only, and the data are from the BaBar collaboration Fzl/ . 

is used instead, the prediction for A{p$) is basically the same. The resummed color-octet 
contributions give A(p^) ~ 1 for nearly the entire p^ range. This is because for our choice 
shape function parameters the color-octet dominates color-singlet over the entire range of 
p. The resummed color-octet calculation agrees with the Belle data as well as the BaBar 
measurement for p^ > 3.5 GeV. However, the BaBar measurement for p^ < 3.5 GeV is 
about 2a below theory. 

Note that the factorization theorem we derive in this paper is only valid in the endpoint 
region and the predictions in the region p^ < 2 GeV should not be taken too seriously. 
The effects that are being resummed are no longer dominant in this region and other color- 
octet production mechanisms may be important. On the other hand, the small value of 
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FIG. 5: The sum of the color-octet and color-singlet contributions are plotted as the upper line. The 
lower line is the color-singlet contribution only, and the data are from the Belle collaboration flty . 

A{p$) could indicate the color-octet shape function we are using is too broad and some 
other mechanism accounts for the observed J/ip production at these lower values of p^. 
Though the results presented here give a fair description of the p^ distribution, the Belle 
collaboration finds that roughly 60% of the cross section comes from the production of a 
J/ijj in association with another cc pair. The theoretical prediction is that this contribution 
is less than 10% of the cross section. Clearly we do not resolve this issue here, and the 
mechanism responsible for the copious production of J/ip + DD must be understood before 
any measurements of the octet shape function parameters can be made at BaBar or Belle. 



18 



2 



-A- 



-Q- 



Oh 
< 












FIG. 6: The angular distribution variable A{p^). The circles represent data from Belle [1L] and 
the triangle data from BaBar PHI /. The solid line is the color-singlet contirbution alone while the 
shaded band includes the resummed color-octet contribution. 



V. CONCLUSION 



Leading order NRQCD calculations of the J/ip spectrum in NRQCD predict an en- 
hancement at maximal energy due to color-octet contributions. Furthermore, the angular 
dependence of the J /if) production due to color-octet contributions is markedly different 
than the color- singlet. Recent experimental investigations H3, H3 of e + e- -> J/ip + X are 
at odds with the leading order predictions. Large Sudakov logarithms appear in the higher 
order corrections to the color-octet production mechanism which invalidate the perturbative 



expansion. In addition, as pointed out by Beneke, Rothstein, and Wise 21], the nonpertur 
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bative v 2 expansion of NRQCD also breaks down at the endpoint, which necessitates the 
introduction of nonperturbative shape functions to resum these large corrections. 

We have studied the color-octet contribution to e + e~ — > J /if) + X using SCET coupled 
to NRQCD. First we showed how the differential rate factorizes into a hard coefficient, a 
jet function containing collinear degrees of freedom, and a shape function. The renormal- 
ization group equations of SCET were used to move all large logarithms out of operators 
and into coefficient functions. This allowed the resummation of leading and next-to-leading 
Sudakov logarithms, and naturally introduced the shape function into the differential spec- 
trum. While previous calculations of Sudakov logarithms in SCET have employed a two- 
step evolution, we did a one-step evolution, similar to the velocity renormalization group in 
NRQCD [3!]. Due to the factorization of collinear and usoft degrees of freedom j25|], these 
two methods give identical results at leading order in the SCET power counting. 

Unfortunately, the shape functions for J /if) production are not known, so our calculation 
is not predictive. In the future it may be possible to extract the shape functions from either 
this spectrum, or from photo- or electroproduction, to be used as input in an analysis of one 
of the other processes. Until that time, the best we can do is to make a phenomenological 
model. We are guided in this effort by knowledge of the scalings of the moments of the 
shape function, which are given by NRQCD power counting rules. 

Using a simple model, we can obtain a qualitative feeling for the effects of the pertur- 
bative and nonperturbative resummations. We found that both the perturbative resum- 
mation or the inclusion of a shape function cause the differential spectrum to be lowered 
and shifted to smaller energies. Together, the spectrum is very broad, extending the color- 
octet contribution to very low energies. Qualitatively, this could improve agreement with 
experimental data for the total rate, the p^ distribution, and the angular dependence of the 
spectrum. However, there are still many puzzles with the J /if) data, in particular the large 
cross section for J j if) produced along with open charm and the exclusive double charmonium 
production |18], both of which exceed existing theoretical estimates. In addition, there are 
discrepancies between the existing measurements by Belle and BaBar regarding the overall 
normalization of the cross section as well as the polarization of the J/if> at largest values of 
p,p. These discrepancies must be resolved and an understanding of the production mehanism 
which gives rise to the large J / ipcc cross section is needed before color-octet shape functions 
can be reliably extracted from the data. 
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APPENDIX A: DERIVING FACTORIZATION THEOREM FOR THE END- 
POINT 

In this Appendix we show that the rate e + e~ — > J /if) + X in the endpoint region can 
be factored into a hard coefficient, a jet function, and a shape function. The jet and shape 
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function involve only collinear and usoft fields, respectively. The analysis is similar to the 
derivation of factorization for b — > X s ^ jiil . 3a |. We begin by writing the differential cross 



section as 



= -^L^^(0|JK0)|J/^ + X)(J/^ + X|J M (0)|0)(27r) 4 5 4 (g-^-^ 



e 



x 

2 



167T 3 S 3 



e 2 



IT / dfye^^(0\4(y)\J/il> + X)(J/il> + X\J^0)\0) 
J X 

- lr , 3^%- (Al) 
where the sum includes integration over the phase space of X. The lepton tensor is 

IT gTlU'to, (A2) 

where pi^ are the momenta of the electron and positron, respectively. 

Matching the QCD current J M onto SCET operators to leading order in A is discussed in 
Section |Tl| Note that Jl{y) picks up an additional phase due to the field redefinition relating 
QCD fields and SCET/NRQCD effective theory fields: 

4(y) = e^nm. y r tM0 B , ± X p + ^ rJg^A • ^ a • o* X - 9 ] . (A3) 

The currents in Eq. ()A3|) are the same that appear in inclusive radiative T decay [26]. To 
lowest order in a s , T^' is the same as in radiative T decay, while rt^/^ is different. 

We will only derive the factorization theorem for the 1 Sq contribution. The generalization 
to the 3 Pj contribution is straightforward. The 1 5'o contribution to the hadronic tensor T^ v 
can be written as 



T / ^4 -i v / i/2(l-x)fi- ?/r (8, 1 So)rt(8, 1 S , o) 

J-iiv — a ye 1 afj, 1 0u 

X 

- r^r^Tg. (A4) 

In the exponent of Eq. ()A4|) . we have used VBJ_ = —y/s(l — r)B^_ so 



q ^Mv^ + P— = V^{l-x)—-^(l-x)—. (A5) 

The term proportional to n M is suppressed by r m 0.08 and can be neglected. We can now 
decouple the usoft gluons in T e g Q using the field redefinition [Hf 

A^ q = YA^ -> W n = YW^Y\ (A6) 

where the first identity implies the second. The collinear fields with the superscript (0) do 
not interact with usoft fields to lowest order in A. After this field redefinition the color-octet 
1 S'o current becomes 

fa*) = V^f^X-P • (A7) 
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The J/ip does not contain any collinear quanta, 2 so we can write 



rp(3a 
J cff 



x i(0|Tr \T B Bf\y)\ ^ |X c )(X c |Tr \T B B^ a (0)] |0> , (A8) 



A',, 



where X u contains only usoft particles and X c contains collinear particles. This is possible 
since interpolating fields for all final state particles can be written with either all collinear 
or all usoft quanta. By introducing an interpolating field, a^, for the J/ip, we can use the 
completeness of states in the usoft and collinear sectors to write 



Eww = L 



(A9) 
(A10) 



Using Eqs. (jSHJ) and (lATol) in Eq. (jSSjl . we get 

1 



rpfia 

1 eS 



d*y e -^/2(i-*)«-»( 1 X l p Y T A Y ty p (y ) a\a^\ Y T A Y t x _ p (0) 1 0) 



x (0|Tr T B Bf )P {y) Tr T^f^O) |0) . (All) 
Next it is useful to define the jet and shape functions. The jet function is defined by 



|Tr 



T B Bf p (y) Tr T U B^ a (0) |0) 



R (0)a, 



Aug 



a/3 



s(l + r) J (2?r) 



e- ik - y J(n-pn-k + P 2 ± ). (A12) 



Note that the jet function is a function of the labels n ■ p and p± of the collinear fields in 
B^ a , and the residual momentum n ■ k. For the process we are interested in n ■ p = —n ■ px 
and p± = 0. Since the jet function is independent of kj_ and n ■ k, the momentum integration 
over these components yields 5(n ■ y) and 5 2 (y±). The identity (valid for q > 0) 



d 4 y e iq - y (0|Tr T B Bf\y) Tr T^B^O) |0) 



21m 



Ae^(0|T(Tr T B Bf )P (y) Tr T^Sf^O) ) |0) 



xB D (0)a, 



(A13) 



shows that the jet function is related to the imaginary part of T— ordered products of the 
composite field B^ a , which is useful for explicit computations. 



2 If the J /ip is produced with large enough energy, the charm quarks could be considered collinear particles 
though the Lagrangian is modified by the presence of the charm quark mass 39j. However, at very large 
energies, for instance at LEP, quark fragmentation dominates over the contributions considered here |40j. 
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The shape function is 



dy- _ 1£+ .(0\xl p YT B Y^ p (y-)ala^lYT B Y^X-^M 



Air 



-e 2 



4m c <<^(%)> 
(OlxtpT 5 ^ a\a^ S(£ + - in ■ D)^ p T B X - P \0) 



(A14) 



The second line of Eq. ()A14|) is most easily derived by evaluating the first expression in 
n ■ A us = gauge (where Y = = 1) , then replacing ordinary derivatives by usoft gauge 
covariant derivatives to restore usoft gauge invariance. The normalization is defined so that 

/ d£ + S (8 ' So) (£ + ) = i-!^-? P V — = 1- ( Al 5) 

This function depends on the state of the cc pair in the J/ip, and thus there will be a different 
shape function for the 3 -Pj 8 ' ) contribution. 

We substitute these expressions into Eq. (jAllj) . then use 



2K, 



1 + cos 2 6) = ^s(l - r 2 )P[r, z\dz . 



(A16) 



where cos 6 is the angle of the J/ip with the e + e~ beam and P[r, z] = a/(1 + r) 2 z 2 — 4r/ (1 — 
r). Note that P[r, 1] = 1. In the phase space we use the short distance mass 2m c not M^. 
The result for the differential cross section is 

^ = a^ 1So) P[r,z] [ d£ + S {8 ' 1So \£ + )j(^(l-r)[^s(l-x + A/M)-£ + }) 



f 1So) P[r, z] [ d£ + S (8 ' 1So) (£ + )J \s(l -r)(l-x + A/M - £ + /M)] . (A17) 



Here we have used the fact that d£ + S^ ,ls °\£ + ) is boost invariant, and 



i4 = -2- 



1- - + (l-x)n^ 

x 



(A18) 



so 



A £ 

sil - r)[] - X+ M-M 



(A19) 



where £ + = £ + Xy/s/M^. In the last line we have expanded p 2 x to lowest order in 1 — x and 
A/M. Finally, writing z in terms of x 



sx + M 2 /x 
s + M 2 



1 — r 



:i-x) 



(A20) 
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where again we have expanded to lowest order in 1 — x and A/M, and defining 

1 + r 



A + M 



1 — r 



;i-o, 



(A21) 



we find p\ = s(l + r)(£ — z) . The factorization theorem can then be written as 



dz 



4 8 > 1So) P[r, z] f di S^\OJ{s{l + r)(C - z)), 



1 + r 

M^l-O + A 
1 — r 



J( S (l + r)(£-*)) 
(A22) 



where 5 ,( - 8 ' l5o - ) as a function of £ is defined such that J d£5 ,( - 8 ' 1 ' So - ) (£) = 1. The integration 
limits on £ are easy to understand: p\ > requires £ > z and the upper limit comes from 
the requirement that da/dz vanish at the kinematic limit z — 1. 



APPENDIX B: THE JET FUNCTION 

In this section we discuss the jet function and its renormalization. We show how to 
renormalize the moments of the jet function and evaluate these moments to 0(a s ). We derive 
the renormalization group equation (RGE) for the moments needed for the resummation of 
endpoint logarithms. We will assume that the usoft fields have already been decoupled by 
the field redefinition, Eq. ()A6jl . and will use the superscript (0) to denote bare fields instead. 

We can invert Eq. (|A12|) to obtain 

J(n-pn-k + p\) = r) (Bl) 

4:71 



xlm 



J d A ye* k -y (0|T (Tr [t b B{{y)\ Tr [T B B ±p (0)]) |0> 



To lowest order in a s , we replace the 5" with A^, so up to a numerical factor the jet function 
is given by the discontinuity of the SCET collinear gluon propagator. Therefore, 



J(n-pn- k + p 2 ± ) = — 1 _ ; Im ^ — - ; 2 - (B2) 



s(l + r) r 1 
Im- — ., 

7r n ■ pn ■ k + p ± + te 

s(l + r)5(n ■ pn ■ k + p 2 ± ) 
6{t - z) . 



In the last line we have used p\ = n ■ pn ■ k + p\ = s(l + r)(£ — z). 

To calculate the renormalization group equations for the moments of the jet function, we 
begin with the definition in Eq. (jBl|) . The moments as defined in Eq. (|30|1 do not correspond 
to local SCET operators, so it is easiest to compute the time-ordered product to one-loop 
then take moments of the result. Since the divergence depends on the moment variable, N, 
we need to introduce an N dependent renormalization for the moments of the jet function: 
J N = Zjfjjtf , where jffl is the A^th moment of the bare jet function. The field appearing 
in the definition of Jn is the composite field B±. For computing the 0(a s ) anomalous 
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FIG. 7: Graphs contributing to Jn up to 0(a s ). 



dimension one can simply ignore the Z factors for the fields in the Wilson lines since these 

>(o) 



counterterms will only contribute at 0{a 2 s ). Therefore, 



%iB K ">, and 



J 



s(l+r) 



N 



4tt 



(B3) 



x / dz z Im 
'o 



J d 4 ye lk - y (0|T (Tr [T B ^ 0)/3 (y)] Tr [t b 5^(0)]) |0) 



where Z% is the wavefunction renormalization of the gluon field, 

z 3 = 1 H C 4 rif- . 

3 4jre I 3 j 3 



(B4) 



The 0(a s ) corrections to the time-ordered product are shown in Fig. In addition to the 
one-loop graphs there are three counterterm contributions. Two come from the Z factors 
in Eq. (jB3|) and correspond to the last two diagrams in Fig. [7| There is also a tree level 
insertion of the 0(a s ) wavefunction counterterm for the collinear gluons, shown as a cross 
in the third diagram of the first line of Fig. This contribution cancels against the term 
proportional to Z3 — 1 in the third line of Fig. [7| 

Therefore, the divergences from the one-loop graphs must be canceled by the diagram 
proportonal to Z^ — 1 in the last line of Fig. [7| Evaluation of the graphs and derivation of 
the renormalization group equation is straightforward. We already calculated the tree-level 
diagram in Eq. (|B2J) and the self energy corrections are standard. The only remaining graph 
to calculate is the second diagram in Fig. which has a divergent contribution of 



s(l + r) 



71 



dz z N s A Im 

2,71 



Px 



id 



i5)p 



(B5) 
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Taking the imaginary part, using 



Im 



fx 



i5 



log- 



1 — i8)p 



TV 



s(l + r) 



— 5(1 — z) log 



s(l + r)_ 



we can easily integrate over 2. Setting the divergent part to zero gives 



ol s C a 



1XE Z 



2vre 



UCa I 2r In, g-g 



(B6) 



(B7) 



Calculating the anomalous dimension in the usual way, we get Eq. (|32|) for the renormaliza- 
tion of the jet function. We can also evaluate the one-loop expression for the renormalized 
jet function from the diagrams in Fig. [7| 



Jjy = 1 + — — log 

Z7T 



+ 



at. 



s(l+r)N J 7r V12 
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C A - -n f ) lo. 



fi 2 N 



s(l + r)N 



+ 



(B8) 



The ellipsis represents contributions which are not enhanced by large logarithms. Note that 
the result of explicit evaluation of Jn agrees with the solution of the renormalization group 
equation to 0(a s ). It is clear that the large logs in the moments of the jet function are 
minimized when ft 2 = s(l + r)N /N. 
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